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The Color Glass Condensate is a universal state of matter which can manifest itself in hadronic processes involving 
small-x partons, like DIS and pp, pA and AA collisions at high energy. Observables are given in terms of multi-gluon 
correlators, whose ensemble evolves according to a RG equation, the JIMWLK equation. We focus on recent progress 
towards its solution which lead to quasi-exact, analytic expressions for the multi-gluon correlators at high energy. 
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In ultra-relativistic hadronic collisions the final state is rather complicated in terms of the type, the number and the 
distribution of the produced particles. Studying multi-particle correlations is one of the main tools used to extract the 
dominant physics in the various stages of such a process. For example, significant attention has been recently given 
to di-hadron correlations in light on heavy hadron collisions, like deuteron-gold at RHIC and the forthcoming proton- 
lead at the LHC. In the former, when moving towards the deuteron fragmentation region, it is observed an increasing 
suppression of the azimuthal correlation of the two hadrons when their transverse momenta are a few GeV ITJEJ. This 
is the kinematic regime which encourages the search for parton saturation in the wave-function of the heavy hadron 
and descriptions of the data based on such a physical mechanism already exist [3] or are towards completion |4|. 

Thus we are interested in hA — > h\hjX, with h a dilute projectile hadron, like a proton, and A a dense target. 
At the partonic level and for example in the case of qg inclusive production, the diagrams are shown in Fig. [T] In 
the first a large-jt quark from the proton splits into a quark-gluon pair which interacts via multiple gluon exchanges 
with the small-x part of the target and then is measured in the forward region, while in the second the splitting occurs 
after the interaction. The target is viewed as a large color field jf a = ^ + a a , the Color Glass Condensate (CGC) (see 
e.g. 10) and the interaction of a parton at transverse position x with such a field is described by Wilson lines like 
vl = Pexp [ig f dx~ a a ( x~ , x)f]. The inclusive qA — > qgX cross section is given by a Fourier transform of 0171 



where S X{X , = (1/N C ) tr(Vj, V X2 ) and Qx lX2 x 3 x t = (W c ) tr(V^ V X2 Vl 3 V Xi ), with the Wilson lines in the fundamental 
representation, correspond to the S -matrix for the scattering of a dipole and a quadrupole, respectively, off the target 
color field. It is easy to understand the counting of Wilson lines in Eq. (fill. For instance, the diagram (a) involves one 
in the fundamental and one in the adjoint and the latter can be expressed in terms of two fundamental ones. When 
multiplied with its c.c. it gives rise to a term with six Wilson lines which is the first term in Eq. ([TJ (plus the second 
term, a l/N^ correction). The rapidity Y is determined by the kinematics of the process and in the forward region 
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Figure 1 : Partonic level diagrams for quark-gluon production in high energy pA (or dA) collisons and in the proton fragmentation region. 



reads exp(-T) — [\k\ expi-y^) + \q\ exp(-y a )]/ yfs. The high energy QCD dynamics is contained in the correlators in 
Eq. Q, where the average is taken with the evolved, at rapidity Y, target probability distribution Wy[a]. This quantum 
evolution of the GCC obeys the JIMWLK equation (see e.g. |8|), which in its simplest form reads 

dW Y [a] 1 r lt / ~,~ ~,~ ~,~ sab 5 5 r n r n 3(6) 



dY 16tt 3 



r , . ah 5 5 d(0)y , ~ 

M mz (l + VtV v -ViV z -VtV v ) — — W T [a]=HW r [a] => -)£- = {HO) Y , (2) 



where M uvz - (u - v) 2 /[(« - z) 2 (z - v) 2 ] is the dipole kernel, the tilde stands for denoting the adjoint representation 
and the functional derivatives act on the upper and lower end-points of the Wilson lines \r and V respectively. This 
form holds only when acting on gauge invariant operators and this is the case for those appearing in Eq. ([l}. 

Two popular ways are used in order to calculate the correlators of interest. The one is to reformulate Eq. |2]) into a 
Langevin equation 0, solve directly for the Wilson lines on a lattice and then average over many "events". The other 
is to construct evolution equations for the correlators (0) Y of interest according to the second part in Eq. Q. Due to 
the structure of the real, the last two, terms of the Hamiltonian this method leads to a hierarchy of equations 1 10]; when 
i/ lea i acts on n Wilson lines, it can lead to n + 2 of them. This second method had prevailed over the last decade, since 
one was mostly interested in the 2-point correlator (S XlXl ) Y in terms of which we can write, for example, DIS cross 
sections and single inclusive gluon production in pA collisions. Eq. |2} leads to the dipole equation {a = a s N c /n) 



d(S Xl x 2 ) Y 
dY 



a C - 



(3) 



whose diagrammatic interpretation in terms of an equivalent projectile evolution is shown in Fig.[2]a at large-A^ c . In this 
limit (and neglecting Pomeron loops flTT)) one factorizes (SS) Y - (S) Y (S) Y to get the BK equation |[T0lfT2ll . which is 
closed and its solution is by now well understood semi-analytically and numerically. The saturation momentum Q s is 
obtained as the transition line from the region where (T) Y — 1 - {S) Y is weak to the region where it is close to unity. 

However, as evident in Eq. ([TJ, less inclusive quantities require the knowledge of higher-point correlators. The 
quadrupole equation lfl3l can also be derived using Eq. |2| and looks considerably more complicated than Eq. ([3]) 



J[(M Xl x 2 z + M XiXiZ - M X2XiZ )(S XiZ Q ZX2X3Xi ) Y + (M XlX2Z + M X2X}Z - M XlX3 t)(S ab&avcsxt)! 



)y a 
dY ~ 4^r 

+ (>Mx2X3Z + A^X3X4Z — -^^XnX^Z )W X3 Z Qx \ ZX4 ) y "I" (J^XiX^Z -^^X^X^Z ~ ^^X} X3z) \" ZX4 Q,X\ X2 X3 z)j 

- (M XlX2Z + M X , X4Z - M XiX}Z - M X2XiZ ){S XiX2 S X3X4 ) Y - (M X]XiZ + M X2X3Z - M XlX3Z - M X2X4Z ){§ X , X J XlX4 ) Y 

- (M XIX2Z + M X3X4Z + M XIX4Z + M X2X3Z )(Q XlX2X3X4 ) Y ]. (4) 

All terms have an easy interpretation; for instance the real term M XlX3Z S ^ Q XlZX3X4 corresponds to Fig.|2]b, while the 
virtual one M XlX3Z § XdXl § XlXi to Fig.|2]c. Even at large-A^ e , where {S Q) Y = {S) Y {Q) Y and Eq. Q becomes a closed 
inhomogeneous and linear equation (with (S) Y known from Eq. ([3])), the large number of transverse variables together 
with the non-locality in the transverse space are prohibitive for the possibility of a numerical solution. We stress that 
the quadrupole is a new object and a priori there is no reason to think that it may be written as a product of dipoles. 

Thus one takes the Langevin approach to calculate multi-gluon correlators numerically, as done for some config- 
urations in [ 14 1 . It was found that the numerical data for the quadrupole are very well-described by a Gaussian mean 



D.N. Triantafyllopoulos / Nuclear Physics A 00 (2012) 



3 




Figure 2: (a) Real and virtual emission of a gluon, a qq pair at \mge-N c , from a dipole. (b) and (c) Two representative diagrams for the corresponding 
emissions from a quadrupole. In all diagrams the dashed line denotes the interaction with the target. 



field approximation (MFA); an extrapolation to arbitrary Y of the McLerran-Venugopalan (MV) model, the typical 
initial condition at some Yq, Such a Gaussian probability distribution Wyla] involves a single kernel, hence all high- 
point correlators can be expressed in terms of the 2-point one. In fact it was suggested long time ago that a "random 
phase approximation" to JIMWLK leads to a Gaussian Hamiltonian [ 15], even though no explicit reference was made 
to high-point functions. A Gaussian ansatz was also used in lfl"6l to estimate l/N^ corrections to the BK equation. 

One can prove that a Gaussian approximation is a quasi-exact solution to the JIMWLK equation ifTTl (TB) . At 
first glance, such a statement looks very peculiar since the JIMWLK Hamiltonian is highly non-linear because of 
the Wilson lines in Eq. Q and which arise from the propagation of an emitted gluon in the background target field. 
This non-linearity is indeed present in the dipole equation but not in the quadrupole one which is linear in Q (at 
large-A^ c ). This already suggests that a Gaussian approximation can be a possible solution with all the non-linearities 
absorbed in its kernel or, equivalently, in the 2-point function. At saturation, where the target is dense, real emissions 
are suppressed and the virtual part H^m composed of the first two terms of the Hamiltonian dominates. This piece 
is evidently of Gaussian form, including the second term; those Wilson lines simply transform the "left" functional 
derivatives to "right" ones which act on the lower and upper end-points of the Wilson lines V and V respectively. 
Now we integrate the dipole kernel over z in the region 1/Q S <s \u — z|, |v — z| <sc \u — v|, with the lower limit imposed 
by our approximation and the upper one chosen to give the dominant logarithmic contribution 2 In [(« - v) 2 2,]- So 
far this treatment leads to a Hamiltonian valid only at saturation, but performing the same approximation in Eq. |3]l 
we see that this logarithm is related to the logarithmic derivative of the dipole w.r.t. Y. With such a replacement we 
arrive at the main result 

H = _J_ f din (S UV ) Y l_S 6_ _6 6_\ 

G VC F J„ V dF W hu 6al Sc^da-J- 

This is a Gaussian Hamiltonian with a kernel which has absorbed the non-linearities and it is most easily determined 
from the BK equation since in the Gaussian approximation one has din (S uv ) Y /dY = (2CF/N c )dln(S uv )^ K /dY. 
Eq. |5]) holds at finite-A^ c , and is correct at saturation by construction and in the dilute limit as can be inspected. 

Using Hq one constructs evolution equations for multi-gluon correlators which have the benefit to be local in the 
transverse plane and are, thus, ordinary differentials equations in Y with F-dependent coefficients. Using a "separabil- 
ity" property of the Gaussian kernel in Eq. <|3j one can write the final solution for a correlator as a local function in Y. 
E.g. the quadrupole at laige-N c with an MV model initial condition (different ones can also be accommodated) reads 

/A . ln[(S 12 ) Y (§ 34) Y /(§ 1 3) Y (§24) Y ] , e . /e . In[<5i4> y <5 2 3>y/<5l3> r <524>y] . . 

In [(S 12) Y (S 34) Y /{S 14) Y (S 23)y\ ln i(S U)y(S23)y/\ S n) Y \ s W Y \ 

(with i = xi). Formally this is the expression derived in the MV model [ 1 3 1 and later on generalized at finite-A^ c 
lfl9ll . We notice that the quadrupole above obeys the "mirror symmetry" (2xi*2*3* 4 )r = (2x1*4*3*2)1'' which holds 
beyond the large-A^ c and Gaussian approximations and is a result of symmetry under time-reversal where time stands 
for x~ Ifl8l . This symmetry is preserved by JIMWLK due to the two types, left and right, of functional derivatives 
and implies that the hadron expands symmetrically in the x~ direction in the course of evolution. Within the Gaussian 
approximation only, the quadrupole is also symmetric under the charge conjugation {Q Xi x 2 x,x4)y = (Qjcixj^x^y- 

At finite-A^ c operators mix and one needs to diagonalize a matrix. The quadrupole mixes with two dipoles, while 
the phenomenologically interesting 6-point operator QS in Eq. ([1} mixes with two more operators and the emerging 



D.N. Triantafyllopoulos / 'Nuclear Physics A 00 (2012) 



4 



0.2 - 
■ 
-0.2 - 



MFA N c = 3 
MFA large- N c 
Factorized N c = 3 
{S)y 

JIMWLK 



(1.2 



0.1 



0.6 



l-(S)j 



Figure 3: (56)yforthe line configuration as a function of 1 -{S) Y . Con- 
tinuous red: JIMWLK with running coupling for six rapidity values from 
Y = to 5.18. Continuous blue: complete MFA (Gaussian approxima- 
tion) result for N c = 3. Dashed magenta: MFA at large-A^. Dotted gold: 
factorizing the average of QS and using MFA at N c = 3 for (Q)y. Dot- 
ted dashed green: (S)y, based on just counting Wilson lines. JIMWLK 
curves from the numerical solution in 1141 . MFA curves are analytical 
expressions in terms of (S), with the latter again provided by 1141 . At 
7 = JIMWLK and full MFA coincide because of the MV model initial 
condition. This remains true at any Y in the dilute and in the dense region. 
A tiny difference, which stabilizes soon, occurs in the transition region. 



cubic equation has a unique real solution. Simple expressions are obtained for some configurations, like the 'line' 
with JC3 = jci and x\ = X2 and thus a single distance determining the value of the correlator, and where we find 

M _ (jV c +2)(jV c -l) /6x fEr (N c + l)(N c -2) 

Using the above we can plot (S(,) Y , defined as the first two terms in Eq. ([TJ, and compare with the numerical solution 
to the JIMWLK equation. As shown in Fig.|3]the accuracy for the full curve is excellent at any value of Y. Needless to 
say, analytic expressions for the multi-gluon correlators are of unparalleled significance when, for instance, one wants 
to performs the Fourier transform of Eq. ([T} and obtain the desired cross section. 
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